Multi-Time Multi-Scale Correlation Functions in Hydrodynamic Turbulence 
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High Reynolds numbers Navier-Stokes equations are believed to break self-similarity concerning 
both spatial and temporal properties: correlation functions of different orders exhibit distinct decor- 
relation times and anomalous spatial scaling properties. Here, we present a systematic attempt to 
measure multi-time and multi-scale correlations functions, by using high Reynolds numbers numer- 
ical simulations of fully homogeneous and isotropic turbulent flow. The main idea is to set-up an 
ensemble of probing stations riding the flow, i.e. measuring correlations in a reference frame cen- 
tered on the trajectory of distinct fluid particles (the quasi-Lagrangian reference frame introduced 
by Belinicher & L'vov, Sov. Phys. JETP 66, 303 (1987)). In this way we reduce the large-scale 
sweeping and measure the non-trivial temporal dynamics governing the turbulent energy transfer 
from large to small scales. We present evidences of the existence of dynamic multiscaling : multi-time 
correlation functions are characterized by an infinite set of characteristic times. 



I. INTRODUCTION 



A comprehensive theory of the Eulerian and La- 
grangian statistical properties of turbulence is one of 
the outstanding open problems in classical physics 
Most studied quantities concerns either measurements 
performed at the same time in multiple positions (Eu- 
lerian measurements) [U, H| , or along one or several par- 
ticles moving with the flow (Lagrangian measurements) 
The latter are optimal to study temporal proper- 
ties of the underlying turbulent flows @, HJ but cannot 
simultaneously disentangle also spatial fluctuations, be- 
n g ba sed on single point, as for the case of acceleration 
9|4l2| . or on a evolving set of scales as for two-particles 
13Ul7| and multi-particle dispersions [l8|, On the 

other hand, neither analytical control nor a firm phe- 
nomenological description of fully developed turbulence 
can be obtained without a solid understanding of the re- 
lation between spatio and temporal fluctuations [20l - l29j . 
In order to access unambiguously spatial and temporal 
fluctuations one needs to set the reference scale and to get 
rid of the large-scale sweeping in the same experimental - 
or numerical- set-up [30H33I ] . The idea here is to exploit 
numerical simulations to define a set of probes flowing 
with the wind, moving on a reference frame sticked with 
a representative fluid particle. In such a reference frame, 
we can access velocity fluctuations over different spatial 
resolution, together with their temporal evolution, with- 



* Corresponding author: cnrico.calzavarini@polytcch-lillc.fr 



out being affected by the large-scale sweeping [34[. The 
interest of these measurements is twofold. First, all at- 
tempts to break the theoretical deadlock in turbulence 
have been hindered by the difficulties in closing both spa- 
tial and temporal fluctuations [U H^, [35| (notice that the 
main theoretical breakthrough in turbulent systems have 
been obtained where temporal fluctuations are uncorre- 
cted [HI). Second, small-scale parametrizations used for 
sub-grid turbulent closure call for more and more refined 
phenomenological understanding of spatial and temporal 
fluctuations [371 ]. 

In this article, we show that quasi-Lagrangian measure- 
ment are able to remove the sweeping effect, revealing 
that correlation times in the Eulerian and Lagrangian 
frame scale differently. Moreover, Lagrangian properties 
posses a dynamical multi-scaling, i.e. different correla- 
tion functions dccorrclatcd with different characteristic 
time scales. The locality in space and time of the en- 
ergy cascade is supported by studying the delayed peak 
in multi-time and multi-scale correlations. The main re- 
sult we present is the confirmation of a theory [2l[ where 
a bridge between spatial and temporal intcrmittency was 
made by means of a refinement of the multifractal phe- 
nomenology (originally proposed for the spatial statis- 
tics of Eulerian velocity in 38]). Energy is transferred 
downscale with intermittent temporal fluctuations, and 
an associated infinite hierarchy of decorrelation times. 
Temporal fluctuations become wilder and wilder by de- 
creasing the scale. 

The article is organized as follows. In sec. II we present 
the numerical methods. In sec. Ill we show the result 
for single-scale multi-time correlation functions, for the 
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FIG. 1: Sketch of the n probes "riding the flow", i.e., at 
fixed positions in a reference frame moving sticked with a 
representative fluid particle x(t). 

bridge relation between temporal and spatial scaling ex- 
ponents and for multi-time and multi-scale correlation 
functions. 



II. METHODS 

A. Sweeping and quasi-Lagrangian reference 

The difficulty in studying the temporal correlations in 
turbulence is associated with the sweeping of small scales 
by means of larger ones. In the case of a flow with a 
"large" mean velocity U, fluctuations u' = u U are 
almost passively transported in space via the Galilean 
transformation x'(f) = x + Ut. This property, dubbed 
Taylor frozen flow hypothesis, is commonly used in ex- 
periments to remap the hot-wire probe readings (which 
are done in time) to measurements in space and is sup- 
posed to be valid as long as turbulence levels are small 
(|u'|/|U| <C 1). In order to study the temporal evolution 
of velocity differences at a given scale r it is hence nec- 
essary to get rid of large-scale (infra-red) effects. This 
can be done by "riding the flow", i.e. by sticking the 
origin of the reference frame on the position of a fluid 
parcel moving in the flow (see sketch in fig. [I}. Such ref- 
erence frame, introduced by Belinicher and L'vov [30], is 
called quasi-Lagrangian. Because of this difficulty, mea- 
sures of the quasi-Lagrangian type have been performed 
only numerically at moderate Reynolds and for small- 
scale quantities [4(| or at high Reynolds in turbulence 
models, such as shell models (where sweeping is absent 



B. Numerical methods 

In the present study we measure multi-scale and multi- 
time correlations in a quasi-Lagrangian reference frame 
from fully resolved high-statistics three-dimensional Di- 
rect Numerical Simulations (DNS) of homogeneous and 
isotropic turbulence. We evolve the incompressible 



Navier-Stokes (NS) equations: 

d t u+(u-V)u = -Vp + ^Au + f, V-u = 0, (1) 

in a cubic three-periodic domain via a pseudo-spectral 
algorithm and 2 nd order Adams-Bashforth time march- 
ing. The forcing f, defined as in [42|, acts only on the 
first two shells in Fourier space (|k| < 2) and keeps con- 
stant in time the total (volume averaged) injected power, 
(f • u)y = const. We report data coming from a set of 
simulations with N 3 = 256 3 and 512 3 , corresponding to 
ReA — 140 and 180 respectively (see Table U for rel- 
evant parameters characterizing the flows). The simu- 
lation e.g. at Kc\ ~ 140 has been carried on for 40 
Eulerian turnover times, T = (3/2)u^ ms /e. We also in- 
tegrated numerically N p = 3.2 • 10 4 tracers evolving with 
the local Eulerian velocity field x(t) = u(x(i),t). At 
fixed temporal intervals we evaluate the fluid velocity also 
at x(t) + ri, with i = 0, . . . , M (spatial distances from 
each tracer). The vectors are chosen to be always along 
one fixed direction, f , and are logarithmically spaced in 
the range between zero and half of the box-size (we use 
AI = 20), sec figure [TJ Similar measurements are done 
also at fixed positions uniformly spaced in the fluid do- 
main. These two set of data are denoted respectively as 
quasi-Lagrangian (L) and Eulerian (E). 



C. Notations and Measurements 

Wc focus our attention on the longitudinal increments 
of velocity differences across a displacement r: 

5u r {x, t) = (u(x + r, t) - u(x, t)) ■ r. (2) 

Notice that we have adopted a unifying notation, for 
us x can represent either a fixed point in space x = 
xo or a point following a fluid particle: x = x(i) = 
J tg u(x(i|xo, to), t) dt + x (a trajectory passing from x 
at time to). We distinguish between the two cases by 
the superscript labels: Suf(it,t) or 5u^(x,t). Note that 
(Suf) p , with overbar denoting time average, is the usual 
Eulerian structure function of order p, furthermore by 
means of crgodicity it can be proved that (5u^) p = 
(6ur) p (therefore for such a quantity E, L labels will be 
dropped). We define the generic multi-scale, multi-time 
correlation functions 1211: 



R ' r[L ' E]y ' (5u R (t))« ■ (6u r (t))P-« 

(3) 

where R and r denote separation vector fixed in space 
and with different magnitude. Note that the L, E dis- 
tinction must be kept for the average of the multi-time 
product in the numerator. Given the correlation func- 
tions we can define an Eulerian and Lagrangian (integral) 
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TABLE I: DNS parameters: N is the number of grid points per spatial direction; Sx — 2-k/N and St are the spatial and 
temporal discretization; v is the value of kinematic viscosity; e the mean value of the energy dissipation rate, ttot is the total 
simulation time; N p total number of fluid tracers; M number of probes at fixed distances from tracer particle; n = (y 3, /e) 

and Trj = (i//e) 1/ ' 2 the Kolmogorov dissipative spatial and temporal scales, u rma = f(iijUj) v / '3) the single-component root- 
mean-square velocity, A = (15 v u 2 ms /e) 1//2 the Taylor micro-scale, T = (3/2)u 2 ms /e and L = u rma T large-eddy-turnover 
temporal and spatial scales; Rex = u rma A/ v the Taylor scale based Reynolds number. 



correlation time as follows [2l|: 
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III. RESULTS 
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Lagrangian correlation time as a function of the Eulerian 
one, a procedure similar to the Extended Self Similarity 
(ESS) [43|], does show a good scaling with slope 0.64±0.02 
in the range [20, 200]^, consistent with 2/3 (Fig. [2]). This 
finding again supports the idea that the limited scaling 
in Figure [2] is due to finite Reynolds effects. 



A. Single-scale multi-time correlation 

We begin discussing the special case of a single-scale, 
multi-time correlation, i.e. R = r. Dimensional inertial- 
range scaling, (Su r ) p ~ (er) p / 3 , provides the following 
estimate for the turnover time of inertial eddies of size 

r, T^ p ~ q \r) ~ r/ ((<5u r ) p ) P ~ r 2 / 3 . On the contrary 
the Eulerian correlation time -due to sweeping effect- can 
be estimated by means of the typical velocity difference 
of the largest eddy, which is proportional to the mean 
square root velocity, <5«l ~ u rrns . One has Tg' p (r) ~ 

r j I (Suj_,) p ) ~ r/u rms ~ r. In the r — > r\ limit both 
correlation times tend to the dissipative scale t v . In Fig- 
ure (inset) we show the behavior of C^r^ir) for both 
Eulerian and quasi-Lagrangian velocity differences and 
for separation scales r £ [2.4, 245]^. On the abscissa 
the time increment r is made dimensionlcss through the 
Eulerian large eddy turnover time T (see Tab HI . We 
clearly see that after a time ~ T all the correlations have 
decreased at least of a factor 50, supporting the quality 
and convergence of our measurements. The main panel of 
Figure [2] shows the integral correlation times both for the 
Eulerian and quasi-Lagrangian case as computed from 
(|4]), in a time integration window [0,T]. The behavior 
is in qualitative agreement with the expected scaling, 
the Lagrangian case being less steep than the Eulerian 
one, however pure power-law scaling seems to be hin- 
dered by finite Reynolds number and system finite-size 
effects. To demonstrate this, we introduce a parametriza- 
tion for the second order spatial velocity structure func- 
tions, with dissipative and large-scale cut-off (see also 
ED): T n (r) = Cl (l + (r/c 2 ) 2 )"/ 2 (l + (r/c 3 ) 2 )-"/ 2 with 
n = 2/3 and n = 1 respectively for the Lagrangian and 
Eulerian case. The parameters c\ , c-i , C3 represent the dis- 
sipative correlation time scale, the dissipative and large 
cut-off scales, respectively. The good quality of the fit, 
shown in Fig. O supports our hypothesis. Plotting the 



B. Intermittency and test of the bridge relations 

It is well known that Eulerian statistics shows inter- 
mittent corrections to dimensional scaling. For example, 
for structure function we have (Su T ) p ~ r^-* where £(p) 
is a nonlinear convex function of p [l[. In 1997, L'vov, 
Podivilov & Procaccia [2lJ provided a possible framework 
to encompass the phenomenology associated to intermit- 
tency also to temporal fluctuations. The idea consists in 
noticing that for time correlations the structure of the ad- 
vection term of the Navicr- Stokes equations suggests the 

relation: T^' p " 9) (r) - rj (jS^f / (Surf 1 ) - r 2 ^ 

[2ll ]. Using the scaling for the Eulerian quantities, 

(Su r ) p ~ r^ p ) one gets to the so-called bridge relations 
(BR) connecting spatial and temporal properties: 

z{p) = 1 - C(p) + C(P - !)■ 

Similar idea have also been successfully applied to con- 
nect the statistics of acceleration and velocity gradients 
[22^ . Plugging the empirical values [H, @ for the Eu- 
lerian exponents in the previous expression, one pre- 
dicts z(p) = 0.67, 0.74, 0.78, 0.80(±0.01) for the orders 
p = 2,4,6,8 respectively. In Figure [3] (main top panel) 
the different moments of the Lagrangian integral times, 
i.e. T^' p q \r) (with p — q = 1), are shown versus the 
scale r. A steepening of the scaling properties with in- 
creasing p can be noticed. In order to enhance the quality 
of the measurements we resort to ESS procedure by plot- 
ting a generic T^ ,p ~ q \r) versus T^' x \r) (see inset of 
Figure [3]) in log-scale. It is now possible to define local 
scaling exponents as: 

z r (p) = dlogT^^/dlogT^ir), 

which according to the BR should scale as z(p)/z(2) for 
r in the inertial range. The result of this local scaling 
exponent analysis is shown in Figure [3] (bottom panel) 
for the orders p = 4, 6. Notice that the BR predicts 
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FIG. 2: Top-inset: Correlation functions of single-scale ve- 
locity differences Cf^'^i e]( t ) as a mnc tion of the time delay 
t, in the (quasi-)Lagrangian (solid line) and Eulerian (dashed 
line) frame of reference. Top-main panel: Single-scale inte- 
gral correlation time T^'^lr) (•) and T^^^r) (□) estimated 
from Lagrangian and Eulerian Cvfr ( T ) correlations. The 
straight lines indicate slopes 2/3 and 1. The fit using the 
function T n (r) discussed in the text with n = 2/3 (solid) 
and n = 1 (dashed line) are also shown. Vertical dotted line 
reports the value A ~ 20?7 for which both correlation times 
are ~ t v . The horizontal line show the integral correlation 
time of longitudinal velocity gradients T^ 1,1 ' (V||ii) = 2.15 r v 
along Lagrangian trajectories, which corresponds to the 
small-scale(ultra-violet) limit, lim, — >o Tj 1 ' 1 ' (r). Bottom main 
panel: Extended Self Similarity plot: T£ M) (r) vs. T^ A) (r) 
and its local slope: d log T^ 1 ' 1 '(r)/d log Tg' X \r) (Bottom in- 



set). The inertial range behavior 
panels. 



.2/3 



is reported in both 



the same scaling properties independently of q. In our 
numerics we find slightly different results for q = 1 or 
q = p— 1. The error bars in bottom panel of fig. gives 
a quantitative estimate of the spread between the two 
results. In the inertial range ~ [10, 100] 77 we find some 
deviation from the K41 values z(p)/z(2) = 1, consistent 
with the BR predictions for p = 4, 6. We notice that 
the predicted intermittent corrections are very small and 
error bars large. Higher statistics and/or higher Reynolds 
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FIG. 3: Test of dynamic multiscaling. Top panel: T ( L q,p q) (r) 
vs. r/rj, with p — 2,4,6,8 and q = p — 1. Inset top panel: 
ESS-like plot of T ( L q ' p ' q) (r) vs. T< M) (r) for p = 4, 6, q = p- 1. 
Bottom panel: local scaling exponents, z r (p), for p = 4, 6 
and 3 = p — 1, 1. Central values are computed as the mean 
between q = p — 1 or g = 1 for each p value. Error bars 
gives the dispersion between the two choices. Horizontal lines 
represents from bottom to top respectively the dimensional 
prediction z(p)/z(2) — 1, Vp denoted as i^41, and the BR 
values z(4)/z(2) = 0.74/0.67 and z(6)/z(2) = 0.78/0.67. 



number may help in giving stronger confirmation to these 
evidences. 



C. Multi-scale multi-time correlation 

We now focus on the most general case of multi-scale 
and multi-time correlation functions in the Lagrangian 
frame. In particular, in the correlation function ([3]) 
we vary the large scale R while the small scale r is 
kept fixed r ~ 77. Note that the velocity difference <5mr, 
precedes in time the difference 6u r ~ v . We are therefore 
interested in the time it takes for a velocity fluctuation 
to cascade down from a large eddy (of size R) to the 
smallest one (of size ~ rf). In Figure 2] (top panel) 

we show the correlations Cj^V l (t) with r* = 2.4 77 
as defined from Eqn. (except for the fact that to 
enhance the contrast instead of Su T we used \5u r \ ). 
The presence of a peak in l (t) for each given R, 
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FIG. 4: Top panel: C [ ^l L (r) for fixed r* = 2.4?? at chang- 
ing R £ [r* : 1207;]. The solid line is a guide to the eyes 
connecting the peak for each given R. Inset bottom panel: 
Cut of C ( 1'^ l (t) at increasing R (from bottom to top) as a 
function of r. Symbol (•) marks the position of the maximal 
correlation value (denoted as T^'^(R)) which increases for 
increasing values of R. Bottom panel: Comparison between 
T ptak( R ) and T i L 1 ' 1) (R)-T ( L 1 ' 1) {r*). The inertial scaling R 2/3 
is also drawn for comparison. 

defines a time, T^^(R), which increases for increasing 
values of R. The presence of the peak can be directly 
associated with the time lapse it takes the energy to go 
down through scales from R to r* , i.e. a direct evidence 
of temporal properties of the Richardson turbulent 
cascade [l|. In the inset of the bottom panel we show 
the curves corresponding to C^'z} l {t), for each different 
R at varying r. 

The scaling behavior of the peak time T^^(R) ~ i? 2 / 3 , 
shown in Fig. 0] (bottom panel), is in agreement with 
what has been measured for Fourier-space based quanti- 
ties by Wan et al. in f34| . 

Also in figureH] (bottom panel) a comparison of T^'^(R) 
with T { L hl) {R) - T { L hl) (r*) (as computed for Fig. 0) 



is shown. It is remarkable to note that the amplitude 
and scaling of T^'^(R) -coming from the multi-scale 
correlation function- is close and compatible with 
T^ a \r) - T^ 1 ' 1 ^?-*) -coming from the single-scale 
correlation functions. This finding provides a clean 
confirmation that energy is transferred down-scale in the 
quasi-Lagrangian reference frame. 



IV. CONCLUSIONS 

We presented an investigation of multi-scale and multi- 
time velocity correlations in hydrodynamic turbulence 
in Eulcrian and quasi-Lagrangian reference frame. Our 
main results are the following: i) We have demonstrated 
that quasi-Lagrangian measurement are able to remove 
the sweeping effect. The integral correlation times in 
the Eulerian and Lagrangian frame are shown to scale 
differently, ii) Lagrangian properties posses a dynam- 
ical multi-scaling, i.e. different correlation functions 
decorrelated with different characteristic time scales, iii) 
Bridge relations connecting single-time multi-scale expo- 
nents with multi-time single-scale exponents are valid, 
within numerical accuracy, iv) The locality in space and 
time of the energy cascade is supported by studying the 
delayed peak in multi-time and multi-scale correlations. 
Temporal fluctuations becomes larger and larger by going 
to smaller and smaller scales, a phenomenon that may 
even affect numerical stability criteria for time march- 
ing, similarly to an effect concer ning spatial resolution 
induced by spatial intermittency |39(. Some issues sim- 
ilar to the ones here discussed have also been addressed 
in a recent numerical study [34| where evidences of the 
Lagrangian nature of the turbulent energy cascade have 
been demonstrated by studying the correlation between 
energy dissipation and local energy fluxes in the quasi- 
Lagrangian frame. While the method followed in |34| re- 
quires a knowledge of the three-dimensional velocity field, 
the approach proposed in the present manuscript needs 
only the knowledge of the velocity at just a few points 
along a Lagrangian trajectory: a measurements which 
may be already accessible in current particle tracking ex- 
perimental set-ups. 
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